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t — . Abstract 
00 

I calculate the leading logarithmic contributions up to two-loop order of the octet 
part of the K° — > tttt amplitude. This sector of the weak chiral Lagrangian is be- 
lieved to be the main source of the enhancement of the 1 = relative to the 1 = 2 
K° — > tttt amplitude, the so-called AI = 1/2 rule. I discuss the procedure of chiral 
extrapolations of lattice data specific to K — > tttt decays and study the implication 
Qh| of the present calculation on these numerically. The latter reinforces the fact that 

Q_|' one has to expect a large enhancement of the 1 = part of the amplitude due to 

^ re-scattering effects between the three mesons. 



■ 1 Introduction 



Chiral logs are introduced during the process of renormalization pp; In the framework of 
dimensional regularization, one has to introduce an energy scale v to ensure the correct 
dimension of observables: Terms like the one on the LHS of Eq. (jl.ljl . generated by loop 
integrals in d = 4 — e dimensions, are re-expressed in the v independent form on the 
RHS 1 . It is necessary to adopt the energy scale v to provide the proper dimension of the 
divergent and finite pieces in a well-defined way. Physically, the emerging chiral logs can 
be associated with the infrared singularities when the masses of the theory approach zero, 



(„-( i„( ))) , (i.i) 



and can produce sizeable contributions to amplitudes. 

Chiral logarithms are an important ingredient in a lattice determination of the K° — > 
tttt amplitude: the lattice community commonly uses a standard method to relate the lat- 
ter to the K° 7T° and K° — > |0) matrix element j2J. However, this approach cannot 

^These terms are multiplied by geometrical factors and polynomials of integer powers in the masses. 
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account for QCD interactions between the mesons, in particular the so called final state 
interactions between the two pions. These effects, though, are expected (at least partially) 
to explain the enhancement of the amplitude with the two pions in the 1 = state relative 
to the 1 = 2 state (A = 1/2 rule) as well as the disagreement of the value of e'/e between 
theory and measurement. CHPT provides a tool to estimate the effects caused by these 
interactions, and the chiral logs are a numerically important part thereof. 
In two flavor CHPT , these contributions are commonly the dominant part of the NNLO 
corrections, which is illustrated in the case of irn scattering in [3]. In the case of three 
flavor CHPT, the double log contributions are however not as prominent. Tabled shows 
the chiral corrections up to NNLO to the pion and Kaon decay constants and the vector 
form factor of Ki 3 jH El E]- The double logs amount to 20 — 35% of the NNLO order 
result, corresponding to around 10% of the total corrections at NNLO to the leading or- 
der result. Since the latter are of importance in the procedure of the chiral extrapolation 
discussed later, the double logs can make a sizeable contribution for such applications. In 

LP NLO NNLO Double Logs 

F n /F 1 0.068 -0.172 -0.050 

F K /F n 1 0.216 0.035 0.06 

f+(0)[K a ] 1 -0.023 0.015 0.004 

Table 1: The chiral corrections up to NNLO for the pion and Kaon decay constants and 
the vector form factor of Ki 3 . The values can however vary considerably, depending on the 
LEC's one employs. The numbers above are calculated with some standard values of the 
NLO LEC's and all the renormalized NNLO LEC's set to zero at y, = 770MeV. 

the framework of CHPT, the interactions between the mesons are encoded in the higher 
order contributions of the amplitude. So far these contributions can't be determined ex- 
actly, since the low energy constants (LEC's) of the weak chiral Lagrangians at Ti-order 
1 (NLO) are not known, let alone the NNLO LEC's. The determination of the leading 
log contributions is thus the only possibility to get an estimate of the corrections one has 
to expect at NLO and NNLO. In [7j a method has been proposed to provide the LEC's 
needed for the K — > nn amplitude from lattice simulations. 
The outline of the paper is as follows: 

In section El I introduce some notation and the chiral Lagrangians used in the calculation. 
Section 01 discusses the procedure of the chiral extrapolation used in lattice determinations 
of the K° — > irir amplitude: The quantity which can be extracted from a lattice simu- 
lation corresponds to the chiral limit of the amplitude, or, more accurately, the reduced 
counterpart thereof. In a further step, the obtained chiral limit value has to be trans- 
formed back to the physical world with non-vanishing quark masses, which in the context 
of CHPT corresponds to the re-introduction of interactions between the mesons (i.e. final 
state interactions), whose importance has already been pointed out. 

This step and the numerics thereof are discussed in section 0] at NLO accuracy, emphasiz- 
ing the role of the chiral logs involved. In section El the procedure is extended to include 
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the double log contributions of this extrapolation, being part of the NNLO corrections to 
the NLO results of section HJ 

In the appendices I outline the more technical aspects of the paper: 

Appendix [X] provides a thorough discussion of the NLO contributions and its approxi- 
mation using only the chiral logs. Appendix |B] presents the procedure used to compute 
the double log contributions of the K — > nn amplitude. The double logs are split into 
different classes of contribution: The "genuine" double logs originating from one particle 
irreducible (1PI) two loop topologies, which can be calculated with the NNLO counterterm 
Lagrangian provided in |Sj. For a second class, the one particle reducible topologies, one 
can employ separate one loop calculations for the respective 1PI subgraphs. The results of 
the needed one loop expressions are given in appendix The last contribution originates 
from lower order diagrams where the masses and decay constants are shifted to their phys- 
ical value, as well as corrections coming from the wave-function renormalization and LSZ 
procedure, up to the required chiral order. The mass and decay shifts are provided in ap- 
pendixlHI Appendix|E]gives a very brief summary of the loop integrals used in appendix O 



2 CHPT Lagrangian 

Here I provide a very brief introduction of the Lagrangians which were used in the calcula- 
tion. For more details about CHPT and in particular the definitions of the building blocks 
used throughout the paper I refer to the numerous review articles, for instance |H| ll()j. 
The lowest order chiral Lagrangian which allows for AS = 1 strangeness changing interac- 
tions is given by: 

£(0) = £(0) + £(0) ^ (2 2) 

with the strong interaction Lagrangian: 

W = + <X+)) > ( 2 -3) 

and the AS = 1 Lagrangian : 

£<J> := CF*g 8 (Au^) ; A := u\ 6 u^ , (2.4) 

where I introduce only the operator which dominates the contribution to the 1 = part of 
the amplitude. &w transforms like (8, 1) under SU(3)l <8> SU(3)r. 

In addition to the lowest order Lagrangians and , I will also use the NLO La- 
grangian HU 112] : 

^-CF^N^W^ , (2.5) 
i=i 

and the NNLO Lagrangian: 

N 

^ = Cg^^wP , (2.6) 
i=i 
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with the respective coefficients: 

NP = (vc)-^Nl 1)r (^e) + ZlA) , (2.7) 

Af } = <jJcr*(^'(jM,e) + 2?A>+(2i? + 2&M)h) , (2.8) 

where Zn(jj) is a coefficient originating from a vertex of t' 1 '. The notation is: 

3 Chiral extrapolation 

For a lattice calculation of the K — > 7T7r amplitude, one normally uses a tree level PCAC 
relation to relate the former to the K° — > tt° and if — > |0) amplitude |2j. The starting 
point is the lowest order Lagrangian, Eq. 12.41 which implies the relation 2 : 

A = A ied . + 0(h) , (3.1) 

with: 

A := (n\°\X^\K°) , 

A,, := -i^^l^l^+^^^l^l/O • (3.2) 
r m z v „ ' 2 .r z A ^ ^ ' 



I use the notation: 



:: ^-red. —'"^led. 



l K- m l ! m 2 :=p K p n , (3.3) 



for the meson mass difference and intermediate meson masses which are evaluated on the 
lattice. A ph corresponds to the physical value of A. 

Since this relation is derived from the tree level Lagrangian, it does not take into account 
any loop contributions, in particular no re-scattering effects between the two pions (fi- 
nal state interactions). However, these physical phenomena are reproduced in a lattice 
simulation, and consequently one has to find a way to split a lattice result into a part 
corresponding to the tree level and the higher order contributions. 

To do so, we write Eq. (J3.1j) again, noting that it is an exact identity if one identifies all 
the matrix elements with the corresponding lowest order quantities (X = Yl^Lo^ 1 ^^)- 

A^ = A^L , (3.4) 
and re-introduces the higher order corrections with appropriate quotients: 

A := qA™ ; A] cd , = qi j£® ; A 2 ^ = q 2 A^ . (3.5) 



2 The A's here shouldn't be confused with the chiral operator used in Eq. (|2.4|l . 
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With these definition one can write ( (A*) = (0|A M |X) = iFp* ; F = q F F 3 ): 

A f „Ql A 1 , lql A\ \ . , 



The expression in brackets on the RHS of Eq. ()3.6|) corresponds to the LEC g$ and is thus 
independent of quark masses. Consequently one is free to take the chiral limit m q — > of 
this quantity, which takes qi, q2 and qp to their chiral values of 1: 

One evaluates the RHS of Eq. ()3.6|) for various quark masses on the lattice and performs 
the extrapolation to the chiral limit, i.e. Eq. (|3.7j) . Examples of how this procedure works 
in practice can be found in [TBI IT3] . 
The rest of this paper will be concerned about q: 

A °° 

q = W) = 1 + £ ^ {n) , (3-8) 

n=l 

which is needed on the RHS of Eq. ()3.7j) in order to bring the chiral limit value of the 
K° — > 7T7T amplitude (i.e. A ph times the expression in brackets) back to the real world 
with the physical (nonzero) quark masses. 

In section 0] I discuss the NLO corrections to q, emphasizing the contribution of the chiral 
logs. Subsequently, in section I will include the double logs, which are a part of 8q^ 2 \ 
into the discussion. The calculation of these is outlined in appendix iBl 



4 The chiral logs of nn at NLO 

Let me start by providing the NLO logarithmic corrections to the K° — > nn amplitude to 
the lowest chiral order result. These corrections were first calculated by Bijnens [To] . 
The K° — > txti amplitude is expanded in its h order: 

oo 

A = J2 HnAin) ■ 

n=0 

At leading order, we have: 

A® = -V2iCF g 8 (m 2 K - ml) . (4.2) 

Altogether six diagrams ( see Fig. 0) contribute to the full NLO result, which is provided 
in Eq. (jC.l|) . The ^-dependence of A^ is drawn as a full line in Fig. ^ 

3 I do not differentiate between and Fk since they will only be relevant in their chiral limit where 
they coincide. 
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Restricting the NLO contributions to the chiral logs using the approximations as outlined 
in appendix El one gets the expression: 

A$ = -V2iCF g 8 ml(l-ml{L w (u) + h K (u)) + ^ m lL n (u) + .^ , (4.3) 

where I used the notation [y: renormalization scale): 

M") := J^pf^ ) 5 X = n,K, V . (4.4) 

The approximation corresponding to Eq. ()4.3|) is drawn as a dashed gray line in Fig. ^ 
closest to the full result. The nearby lightly dashed line corresponds to the leading order 
result plus the full logarithmic NLO contribution. 

q 
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Figure 1: The renormalization scale dependence of q = \A^ \/A^ in terms of Fq. The 
full line corresponds to the full expression, Eq. dHU), with all LEC's N^ r (v) set to zero. 
The lightly dashed gray line is the approximation taking at NLO only the (full) chiral logs 
into account. The dashed gray line below, closest to the full result, corresponds to the 
approximation given in Eq. (|4.3j) . 



In order to get an estimate for = \A^\/A^ , I will use the natural scale of CHPT, 
47r F^l \fWf ~ m p , Nf being the number of dynamical fermions where one expects 
the values of the LEC's not to have acquired too large values due to running. Using F as 
input, one gets: 

gW(i/ = 0.77GeV) ~ 1.52 (1.61) , (4.5) 

which corresponds to the full amplitude, as given in appendix0 with all NLO LEC's set to 
zero; The number in brackets is the value of the approximation given in Eq. (J4.3|) . Please 
note ( Fig. Q ) that the latter simple expression is able to mimic the full amplitude ( with 
iVf )r (z/) EE 0) With an accuracy well below 10%. Since the lack of knowledge of the LEC's 
introduces an error of the same order or probably even higher, it does not really make a 
big difference if one uses the simple form of Eq. (J4.3|) instead of the full expression, Eq. 
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5 Chiral extrapolation including the double logs 



In this section I include the double log contributions to A into the discussion of the factor 
q. The more technical details of how to calculate these is outlined in appendix El Using 
this the approximation, we have: 

A = A (0) + &A (1) + SASl , 

= A® (1 + 6® + $) . (5.1) 

I will again work with the approximations for the chiral logs as outlined in appendix El 
i.e. Eq. fTOjl and ifPjl . 

In this scheme, the following terms contribute to the amplitude: 

gj 0") = <(4|^(v)L^(r/) + ag^H 2 + «£UM(v) 2 ) , (5.2) 

where I introduce an intermediate meson tum for the double log terms which were calculated 
via the leading two-loop divergences. For these one cannot determine whether pions or 
kaons in the loops generate the logs. The first two terms in Eq. (|5.2j) originate from 
diagrams where we performed the loop integration explicitly ( Terms proportional to L 2 
are not generated). In the following I will parametrize the intermediate mass M by: 

m M (t) = m\mY t E [0,1] . 

Similarly as for the NLO result, Eq. ()5.2|) is not well defined in the limit where m n — > or 
m M ~^ 0. One should again view Eq. ()5.2j) only to be an allowed approximation for the 
double logs for the physical (nonzero) masses. 

The term m%5^ K L^{y)LK{v) in Eq. ()5.2jl has a similar origin like the analogue piece of the 
NLO result, while the third term, m^a^ > M L M (z/) 2 , is generated by graphs with 1PI two- 
loop subdiagrams. We will see below that the double log contribution is rather sensitive on 
the value chosen for tum- The 1PI two-loop diagrams contributing to this term are either 
of the "eight" or "sunrise" topology. For the former, the loop integrals are factorized, 
leading only to chiral logs of the form m 2 x Lxi from which follows that tum — m K or t = 
for this class of diagrams. For the sunrise topology, however, terms of the form m\ c L 2 K 
and mj^L^Lx can be generated [HI Ej, presumably shifting £ to a nonzero value. I am 
not aware of a method how to accurately estimate the value of t. However, it was already 
noted by Bijnens et al. in 0, discussing the double logs in the strong sector, that a large 
value of t leads to unrealistic large double log contributions ( They used t = 1/2 ). In the 
following I will use the conservative values t G [0,0.1]. The inability to estimate t in an 
exact way unfortunately hinders us in making very accurate predictions for q. 
The double logs originate from four classes of diagrams, as outlined in appendix [HI The 
separate contributions of these diagrams as well as more technical details can be found 
there. Here I provide only the final result: 

/ 1 27 

A ( 2 >) = -V2iCF g 8 m 2 K (l - m 2 K (L^(u) + -L K {v)) + — m*L» 

+rr £( _ ™LJu)L K (y) - -^-L K (is) 2 + ™L M (^) 2 ) + , (5.3) 

K\ 10g n\ ) K\ ) 12g60 A V ) T M\ , ) T J , V ) 
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for the amplitude given in terms of Fq 4 . 

Fig. |21 shows the v dependence of q using the full NLO amplitude (N^(u) = 0) plus the 
double log contribution. Fig. El illustrates the accuracy of the approximation, Eq. (|5.3|) . 
based on the use of Eq. ()A.2j) and ()A.3|) . by comparing it with q calculated with the "full" 
result for the double logs. Fig. |U displays the dependence of qi on the intermediate mass 
m M parametrized by t. 

As already noted in [T3|, the ^-dependence of q is quite strong. This dependence would be 
canceled by the LEC's, whose values are unknown and consequently set to zero at all values 
of v, making exact predictions impossible. Employing the same reasoning as in section^ 
I will use the renormalization scale v = 0.77GeV in order to get an estimate of q. Various 
values of q at v = 0.77GeV at different values of t are provided in Table El 





q{v = 0.77GeV) 


NLO 






1.52 (1.61) 


NLO + Double Lo c 


is (t = 


= 0.00) 


1.64 (1.75) 


NLO + Double Lo c 


r S (t = 


= 0.05) 


1.77 (1.88) 


NLO + Double Lo^ 


5 S (t = 


= 0.10) 


1.91 (2.03) 



Table 2: The "enhancement factor" q at NLO with and without the inclusion of the 
double logs at v — 0.77GeV, using F as input. The numbers correspond to the full NLO 
expression provided in appendix O with Nl(v) = 0. The numbers in brackets correspond 
to the approximation as given in Eq. ([5.3)1 . 



6 Impact on present lattice results 

In this section I discuss the implication of the presented results on the values for the Re(y4 ) 
provided from lattice simulations. 
The experimental value of Re(Ao) is: 

Re(A ) cxp ' = 33.3 ■ 10~ 8 GeV (6.4) 

The CP-PACS collaboration calculated Re(Ao) using quenched domain wall fermions [T3*] . 
In their paper they use the reduced matrix element, Eq. ()3.2j) . and extrapolate it to 
the chiral limit, analoguesly to Eq. (|3.6[) . For this means they either use a quadratic 
polynomial (Fit I: Re(A ) = £o + £,i m \i + ^("^m) 2 ) or chiral logarithms (Fit II: Re(v4 ) = 
£o + + £,2m 2 M L M )- The outcome is: 



Rg(^-o) chiral limit 



16.5 • 10~ 8 GeV Fit I 
20.7-10~ 8 GeV Fit II 



4 Fq should also be used in the definition of the Lx , Eq. I|4.4() . 
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Figure 2: The renormalization scale dependence of q at NLO (full line) and at NLO with 
the double log contributions included with an intermediate mass rriM(t = 0) (lowest dashed 
line), m M (t = 0.05) ( intermediate lightly dashed line) and m M {t = 0.10) (uppermost 
dashed line) respectively. All iV7(z/)'s are set to zero. The most natural choice to get an 
estimate of q is to set all LEC's to zero at the natural scale of CHPT, v ~ m p , which results 
in the values 1.52 (NLO), 1.64 (t = 0), 1.77 (t = 0.05) and 1.91 (t = 0.1) respectively. 



They do, however, not implement higher order CHPT corrections to the K — > tttt amplitude, 
i.e. they use q — 1. If we take q = 1.77, corresponding to v = 0.77GeV and t = 0.05 ( see 
section |S] ) , one gets the values: 



Since the above numbers can only serve as an order of magnitude estimate, I do not pro- 
vide any error bars for them. The errors for q are probably in the range of 20% or even 
higher. The important statement one can however extract from the above numbers is that 
the corrections due to q can explain the discrepancy between the experimental number and 
the lattice value (which corresponds to the chiral limit value of Re(ylo) )• 
The RBC collaboration ^3] which also used quenched domain wall fermions in their sim- 
ulations account for the chiral logs of the K — > 7T7t amplitude, but use the old result of 
Bijnens jT^l, which did not take into account all the diagrams contributing to these and is 
therefore incorrect ( see section El). Their result with and without the "old" chiral logs is: 




Re(A ) 



RBC 



20.9 ■ 10~ 8 GeV chiral limit 

29.6 • 10~ 8 GeV with "old" chiral logs 



If we again implement the factor q = 1.77, we get: 



Re(A ) 



RBC 



37.0GeV , 



(6.5) 



which agrees well with the experimental value within the rather large errors one has to 
assume. 




Figure 3: The v dependence of q for 
m-M^ = 0) = m K using the full 
result of the double log contribution 
(lightly dashed line) and the approxima- 
tion given in Eq. (|5.Hjl (upper dashed 
line). The full line correspond to the 
NLO result. 

7 Conclusion 
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Figure 4: The dependence on q at v = 
0.77GeV on the meson mass value used 
in the logarithm Lm in Eq. (|5.3[) . 
parametrized by tum — rn^m 1 ^ 1 . The 
full line corresponds to the NLO value. 



I have calculated the contributions of the (leading) chiral logs at NLO and at NNLO us- 
ing Renormalization Group methods. I studied the result numerically and confirm and 
even reinforce the known result that these contributions are quite large ^Hj- With the ap- 
proximations and assumptions used in this paper, the lowest order K — > tttt amplitude is 
enhanced by a factor between 1.64 and 1.91 due to loop effects, depending rather strongly 
on the average value employed for the mass of the mesons in the loops generating the logs, 
but in all cases further enhancing the NLO value of 1.52. The numerical results provided 
here thus point in the right direction to explain the AI = 1/2 rule as a low energy QCD 
effect. 

The above given numbers can only serve as a rough estimate of the higher order contribu- 
tions, since we are not in a position to include the contributions of the LEC's at NLO, and 
other additional contributions at NNLO. However, given that it is not possible to calculate 
the higher order LEC's in this sector so far, these estimates are yet the best one can do 
with an acceptable effort to get an idea of the correction to the lowest order result one has 
to expect. 

Even though not providing accurate numbers, the calculation presented here shows that 
one has to anticipate rather large higher order CHPT corrections to the lowest order result 
of the K -> ( 7T7r)/ = o amplitude. This again indicates the importance to make further ef- 
forts to include these contributions, physically corresponding to the effects caused by the 
interaction between the three mesons, in an accurate way. 
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A Detailed discussion of the NLO contributions 

In this section I give a detailed discussion of the NLO corrections and the approximations 
I made. The diagrams which contribute are drawn in Fig. 

Below I will work with the approximation where one restricts the NLO corrections to the 




a b 




d e f g 



Figure 5: Diagrams contributing to the K° — > tttt amplitude at NLO: a: Leading order, 
b: Wave-function renormalization, c: Decay constant renormalization, d: first tadpole, e: 
second tadpole, f: first unitary correction, g: second unitary correction. Not drawn are the 
two counter-term diagrams. 

chiral logs: 

A = A<®(l + &A£ + ...) , 

= A {0) (l + a n L 7T (u) + a K L K (u) + a v L ri {u) + ...) . (A.l) 

The ax above are polynomials quadratic in the meson masses. 

In addition I neglect all m 2 contributions relative to m 2 K in the polynomials (a's) and write 
the resulting quantities with a tilde: 

a x = m 2 K a x + 0{ml) . (A.2) 

This approximation will be particularly useful in the case of the NNLO contributions. 
Furthermore I identify the logarithms: 

L K (u)~L v (u) . (A.3) 

All following quantities for which the above approximations, Eqs. (jA.2|) and (jA.3|) . are 
employed will be written with a tilde. 
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Using these, I get for SA 



= m 2 K (-LM-±L K (v)) . (A.4) 

This result corresponds to the lowest order amplitude in Eq. ([4.2)1 . expressed in terms 
of the bare decay constant 5 . In the following I will write all quantities in terms of F ; 
Expressing amplitudes in terms of the renormalized decay constant F n increases generally 
the v dependence artificially. I use the values: 

F w = 0.0924 GeV ; = 1.06 , (A.5) 

which should be precise up to two loop order [TH]. Using this approach, the results will 
not get distorted by omitting known contributions to F n like LEC's etc. . 
Please note that Eq. (jA.4j) does not have a well defined limit m n — > 0. The above 
approximation is meant to be used with the physical values of the meson masses. If we 
leave this regime and take the m n — > limit, one has to include the associated J function 
in order to cancel the occurring divergence. In principle it would be aesthetically more 
appealing to have an expression at hand that respects the underlying chiral symmetries 
and limits, which could be accomplished by re-introducing parts of J contributions in Eq. 
(|A.4[) . I will however limit ourselves to the strict use of the chiral logs as unambiguously 
defined in Eq. ([1.1)1 . since this approach defines a well-defined setting. 
The separate contributions of the diagrams to a n , olk and a v are given in Table El 
In order to obtain an expression which agrees with the full logarithmic corrections at the 



Diagram 


a 


b 


c 


d 


e 


f 


g 


total 



















-1 





-1 


a K 





7/12 


1/2 


-2 


1/3 


- 1/6 


1/2 


-1/4 







1/6 





-22/27 


8/27 





23/27 


1/2 


Sum 





3/4 


1/2 


-76/27 


17/27 


-7/6 


73/54 


-3/4 



Table 3: The various contributions of the logs of the separate diagrams. 
5 — 10% level, I also include the most dominant subleading contribution 6 : 

A ( ! g } = -V2iCF g 8 ml(l-ml(L n (u) + h K ^)) + ^mlL 7T (u) + ..}) . (A.6) 
In the first paper treating these corrections the following value was provided [15j : 

5 Replacing Fq with F n shifts the constants in front of the chiral logs. 

6 This is necessary since the coefficient of the subleading term is a factor 3 times larger than the sum of 
the leading coefficients. 
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Writing the result, Eq. (jA.4|) . with the above normalization, using L n ~ Lk — Lm, in 
terms of F and F n respectively, one gets: 

\m 2 K L M iyj) , (A.8) 

l j-m 2 K L M {v)) . (A.9) 

Equating L n and Lk I found numerically to be a rather crude approximation; I have 
checked that this identification results in values which disagree by approximately a factor 
of two from the full result of the logarithms. 

In the original paper ^Hj the diagram e in Fig. |H]is missing. This together with a different 
treatment of the m^L^ contributions generated by diagram f is the origin of the above dis- 
crepancy 7 . Let us point out in this context that I checked my full NLO K — > nn amplitude 
on which the results above are based with the output of the Mathematica package provided 
in and found complete agreement. 

B The calculation of the K —> tttt double logs 

In this section I provide the details of the calculation of the double logs of the K — > 
ixix amplitude: 

A = (7r(p 1 )7r(p 2 )(J = 0)\X$ 1 \q w )\K(q K )) . (B.10) 

Throughout the calculation, I allowed the weak Hamiltonian to carry the momentum q w , 
a setting which can be used to account for the final state interactions dispersively [2*U] . 
However, since the resulting expressions get too cumbersome, I will only provide the results 
for the physical case, q w = 8 . As in the previous sections I will exclusively consider the 
octet part of the nonleptonic weak chiral Lagrangian. 

The double logs of ft-order 2 contributing to Eq. (|B.10|) are produced by two-loop diagrams, 
each loop contributing one logarithm. The amplitude can be expanded in the border : 

A = A^ + A w + A {2) + 0(h 3 ) . (B.ll) 

All "genuine" double logs originating from graphs with a two-loop 1PI subgraph can be 
calculated with the 4j poles in the two-loop counterterm, depicted in FigJHl diagrams a and 
b. Additionally, there are graphs with two independent one loop subgraphs, corresponding 
to diagram c. A further contribution originates from LSZ/wavefunction renormalizations 
of graphs of lower chiral order, O(jo 2 ,p 4 ). Finally, the lowest order masses and decay con- 
stants of these lower order graphs have to be shifted to their renormalized physical value, 
up to the required chiral order. This two last contributions are represented by diagram d. 



7 private communication with the author of |15| . 
8 The full expression can be obtained from the author. 
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Figure 6: Diagrams contributing to the K° — > nn amplitude at NNLO. The square boxes 
represent 1PI loop graphs with one insertion of a weak vertex and loop order given inside 
the boxes. The circles are the analogues with no weak vertices, a: genuine 1PI, b: genuine 
1PI subgraph attached to a strong tree vertex, c: A weak and a strong subgraph with loop 
order 1 respectively, d: Lower order graphs (represented by a box with no associated loop 
order) which experience mass/decay constant shifts or LSZ renormalization. I have not 
drawn the associated counterterm diagrams. However, for diagram a and b I calculated 
the double logs via their associated counterterms, and for diagram c and d I used partially 
the same approach. 



For the graphs which consist of two one-loop subgraphs, diagram c, one can calculate 
the subgraphs separately ( One- loop Ktt scattering and K — > |0) ). However, there is a 
subtlety one has to account for: Since the Kaon which propagates on the one-particle- 
reducible line of these graphs is offshell, one cannot cut the graphs and identify the two 
resulting subgraphs as physical processes, since for this to be true, the Kaon would have 
to be onshell (an asymptotic state). As a consequence, one cannot treat the "genuine" 
graphs a and b, the 1PR graph c, and the LSZ/wavefunction contributions, diagram d, 
separately. In particular, only the sum of these graphs will cancel all divergences. 
For the graph c, I calculate the double logs explicitly by multiplying the two separate one- 
loop graph expressions, and get a residual 4? divergence. This divergence, however, can 
be absorbed into the counterterm of the remaining, j% generating graphs. If we proceed in 
this way, the shifted counterterm will cancel the 4j divergences of the rest of the diagrams. 
For those one can calculate the double logs as follows: All loop divergences have the 
structure ( see Eq. (jl.ljl ): 

Q = z/- £ (A-^/(^) + ...) , (B.12) 
which ensures scale independence. I used the definition: 

(2) 

If we define a basis W$ of operators for the border 2 Lagrangian, the sum of the diver- 
gences of all "genuine" two-loop diagrams will result in the structure: 

N N 

Cg 8 J2 a i W i 2)Q2 = Cg^- 2e Y,^ W i 2) ( A * "2A/iAfM + Vm{v? + -) .(B.13) 
i=i i=i 
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Note that one has generated a nonlocal divergence ALm- These structures will, however, 
be canceled by the contribution of one loop diagrams with one insertion of a vertex from 
H^ 1 ' . Thus, the sum of the divergences of all diagrams will be free of nonlocalities and be 
canceled by the Lagrangian £/ 2 \ Eq. ()2.6|) . 

Since the double logs of the loops always show up in the combination A 2 + fi 2 M , see Eq. 
()B.13|) . and the A 2 divergences have to be canceled by the Zf 2 A 2 term in Eq. (|2.8|) . the 
double logarithms can be calculated as follows ( The £ 22, s correspond to the 1/e 2 poles of 
the sum of all loop contribution, £ 22 = —Z 22 ): 

= -Cg % v-^A\ m {vf . (B.14) 

Note that of the chiral order p 2 and p 4 graphs, which are affected by wavefunction and 
LSZ renormalization, only the tree graph will generate ^ terms: the one-loop graph is 
finite and the wavefunction and LSZ renormalization shifts on the one-loop graphs are of 
border 1 and can therefore only introduce ~ divergences, and consequently we can discard 
the one-loop graph from the list of terms which contribute to \ divergences. 
The final result of this calculation, split into a sum of contributions originating from the 
different diagrams drawn in Fig. El is provided in Table EJ 



Diagram 


~(2) 


~(2) 
a KK 


~(2) 
a MM 


a 








18599/1944 


b 








-1213/972 


c 


-17/27 


269/540 


4/3 


d 


- 17/12 


-1561/2592 


181/81 


Sum 


- 221/108 


- 1349/12960 


7703/648 



Table 4: The various contributions of the double logs of the separate diagrams, which are 
drawn in Fig El 



C One loop amplitudes 

For completeness I provide here the full one loop K° — > -kit amplitude as well as the 
K°K° — > 7T7T and K° — > |0) amplitudes, needed in intermediate steps of the double log 
computation. All these amplitudes have been calculated with the generalized kinematics, 
where the weak Hamiltonian is allowed to carry momentum. The results displayed below 
correspond however to the amplitudes with the weak Hamiltonian at rest. The general 
result can be obtained from the authors. 

The functions used in the amplitudes given below are defined in appendix lEl 
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The renormalization scale dependence of the LEC's Ll(v), N[(v), which cancels the in- 
dependence of the chiral logarithms ^x{v) will be suppressed in the following. All calcu- 
lations have been performed with FORM |21j . 



C.l K° 7T7T 



(ttVI^ 1 )^ ) 

7!^ ( - 2 + jk< 9< + 2 <) + (-T + ^ + »* + (i3 - 2 3^ 



iCF v9B A r n , N, 8 2 , n _ 2 , , , 27 , , , , ,37 m 2 ^ 

r? ""77 

i 4 2 
777 777 

-jflAr " A,) + 4-^f(8^ + 6LI - 2NI - 4iV? - JVJ + 2iV[ + 4^ + 2iV[ 2 ) 
' 4< (16LS + 2Z7 5 - Nl + 2JVf - 2NI - JVJ) - ^^ ^(O - ^^«) 



, 2 4 , — 5 ^.5 , — , 7 ^.yy 2 7T7T \ j\ / QfP2 

7T 7T 7T 

+ 2mH) + ^m A K JK,(ml))) + 0(^ 2 ) • (C.l) 
I checked the above result with the Mathematica program recently provided by ^H] • 



C.2 K°K° 7T7T 
(ttttI^IK ^)/? ^))^ 



1/1 2 , a / 13 2 17 2 \ 2 i ™ ,1 2 2 ' 4 \ 

■ 7^ m K + A (— r< - — m K ) m x + — (-m w m^ - — ) 



2F 2 V3 v 90 " 45 -f 12 9 w 12 

- m 2 ^ ml- ,2 lm 2 K ., 2 ,11 lmi 1 m% N 2 

^ - + "X<3 - S^" 4 + "* ( 30 - 6^ + 6 

/ 1 o ^ 2 ^^7T / ~^ 2 ~^ 2 \ \ r r / 2 2 8 4 \ 20 „ 9 2 

+32L r 6 mlm 2 K + L^-m* + 32m 2 n m 2 K + |mjj-) + J^(m^)(^m 2 r - ^m 2 K )m 2 K 
C.3 K° -> |0> 

(0|^ 1) |^°(gx))^ =0 = 

lg ^ 8 ( - y Am^A + - ^ K m 2 K + ^ (m 2 - 477-4) 

-8A( + 2A[ m 2 f + N r n (2m 2 K - m 2 ) + 8A 2 >4)) + 0(h 2 ) . (C.3) 
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Please note that the result given above correspond to the K°K° — > nir and K° — > |0) amplitudes 
evaluated at the unphysical point where the momentum of one of the Kaons vanishes. I 
checked that I agree with the results of [22] for the K°K° — > 7r7r if I use the physical 
kinematics. 



D Renormalized masses and wavefunction renormal- 
ization 

In subsection ID. II I provide a list of the logarithmic contributions of the shifts of the 
bare masses to their renormalized value, up to NLO. In subsection ID.21 the logarithmic 
contributions of the wavefunction renormalization shifts are given up to two-loop order. 



D.l Masses 

For the masses I only need the h = 1 corrections, since the tree Lagrangian given in Eq. 
12.41 generates only dynamical masses: 

m \ = mf 2 (l + 5m 2 x ) ; X = ir,K,r]. (D.l) 

The log contributions of the NLO shifts thereof are [2*3] : 



ml /n 2 ml. ,7 ml l6m 2 K 
m^ 6 m* 9 m^ 9 mi 



D.2 Wavefunction renormalization 

Z X = 1 + 8Z X ; X = Ti,K,r l 

with: 

oo 

5Z x = ^h n 5Zf . 



n=l 
(2) 

We split SZ^' los into a piece linear (5i) and quadratic (62 ) in the logarithms: 



The one-loop contribution is 



A7 W 4 „ , 2 „ 



S7 (l) 1 . 1 

6z k\ os = 2^ + /ix+ 2^ 
8Z^ oz = 2 fix ■ 
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From SZjp one only needs ^^xL f° r ^ ne calculation. Since the wavefunction renormaliza- 
tion at this order is not available in the literature so far , I use the counterterm Lagrangian 
provided in [21] to obtain it. As already discussed in appendix El one has to assume an 
unspecified intermediate meson mass tum m the loops generating the chiral logs: 

(2 ) _ /113 m* 53 m|m| 73m^ , 
r 7(2) _ ^17 43 m^m^ 32 , 2 

7 {2) /5 \m\m\ 19 m 4 K 



x y( 2 ) — ""tt x " l -k'"'K 1 x ^ ""jf \ ,,2 



^m M 6 77% 6 ?r% 



The above expressions are not well-defined in the limit niu — > 0, and one should only 
employ them for physical masses, i.e. m M > (See section HJ). 



E Loop integrals 

All tadpole diagrams translate into the integral: 



d d l 

(2n) d P - m\ 
1 „^ d s 



( 47r )d/: 



r r(i - -)m~ 



1 , ,m? 



*(-2A + ^ln(-f)) + 0( £ 



where I use: 



e:=4-d ; A := Ne , (E.2) 

z/ being the renormalization scale and c parameterizing the renormalization prescription 
( c = — (log(47r) — 7 + l)/2 in MS). The simplest integral which occurs in the diagrams 
corresponding to the unitarity corrections is: 

Jab{q2) := _1 / (27p ¥^ml (I - q)" - m\ • (K3) 

Similar integrals to J with a polynomial in / in the numerator can be written in terms of 
a combination of A and J. J can be split into a divergent and finite part as follows: 

J ab (t) = ~(2A + 2k ab ) + J ab (t) , (E.4) 

with 



_ N \og(m 2 a ) - \og(m 2 b ) 

^ab n 9 2 

2 mi — mi 

a b 



(E.5) 
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and 



■hbit) 



Jab(t) ~ Jab(0) , 



(E.6) 





2t 



A 



log( 



(t + A) 2 - A 2 
(t - A) 2 - A 2 



) 



; A < . 




with: 



E : 




A : 



\/[A|; X: = (t-(m b 



m a) 2 ) (* - (m a + m h ) 2 ) 
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